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Abstract 

For the case of finite life-time broadening the standard Kubo-formula for the optical 
conductivity tensor is rederived in terms of Green's functions by using contour integra- 
tions, whereby finite temperatures are accounted for by using the Fermi-Dirac distribution 
function. For zero life-time broadening, the present formalism is related to expressions 
well-known in the literature. Numerical aspects of how to calculate the corresponding 
contour integrals are also outlined. 



1 Introduction 

A general theory of linear response functions has been formulated several decades ago @, [[| 
that became the basis of numerous studies of correlation functions, susceptibilities and transport 
properties in solid matter. Although early attempts were already devoted within this theoretical 
framework to magneto-optical phenomena, i.e., the Faraday and the Kerr effect, it was the 

rapid development of both, band structure methods and computational facilities that brought 
about the renaissance in this field from the end of the 1980's (see ||, |3j and refs. therein). This 
progress was built on the assumption that the original expression for linear response functions, 
derived on the basis of a non-relativistic approach, is still valid when applying relativistic band 
structure methods. Recently, this has indeed been rigorously confirmed by Huhne and Ebert 

i- 

Within an independent particle picture the interband part of the optical conductivity tensor at 
finite temperatures, £ M „(u;) is given by [|l|], 



sn f \ h \ " /( e n)(l /( e m)) / ^nm^mn ^nm^mi 



iV r^n En "Cm {hu + € n - € m + i5 KuJ + 6 m - € n + id 

where u is the frequency, /(e) is the Fermi function, V is the volume of the system, the 
Jnm = ( n I I m ) are t ne matrix elements of the current density operator (// = x, y, z) 
with respect to the eigenf unctions | n) of the one-electron Hamiltonian corresponding to energy 
eigenvalues e n . Life-time broadening effects due to different scattering processes are partially 



accounted for by a finite value of 5 = h/r ~ 0.03 Ry, where r can be associated with a 
relaxation time 0, [|. Such a life-time effect can naturally be incorporated into linear response 
theory either by assuming an adiabatic switch on of the external field 0(r, t) = 0(r)e t//T || or, 
equivalently, by convoluting the corresponding frequency dependent conductivity for 5 — > +0 
by a Lorentzian of a finite halfwidth 5. 

Expression ([I]) is usually directly evaluated by applying standard band structure methods for 
ordered bulk (three-dimensional translational invariant) systems by associating the states | n) 
with (three-dimensional) Bloch-functions 0,0. This kind of approach, however, cannot be 
extended easily to layered systems or systems with a surface, i.e., to systems with (at best) two- 
dimensional translational symmetry. Furthermore, it is totally unsuited to deal with disorder 
such as inter diffusion effects in magnetic multilayer systems. 

For zero life-time broadening an alternative approach consists of splitting Eq. ([I]) into an ab- 
sorptive and a dispersive part, which in turn are interrelated via the Kramers- Kronig relations 
0,0. Thus it is sufficient to directly calculate the absorptive part only, the corresponding 
Kubo-formula of which can be expressed in terms of Green's functions as was shown by Butler 



for to = 0, and was used recently by Banhart |10| for to > 0. Even much earlier, Luttinger 
H derived a general formula of the conductivity tensor which builts on the Green's function 
of the system. However, when using this method the corresponding energy integrals have to 
be performed along the real axis which makes a calculation of S M ^(o;) enormously tedious and 
also numerically hard to control. 

In the present work a suitable extension of Luttinger's formula is given in terms of contour 
integrations that avoids the above mentioned numerical difficulties, is suitable for applications 
to magnetic multilayer systems and for treating alloying effects (interdiffusion etc.) by means 
of the Coherent Potential Approximation (CPA). 



2 Contour integrations 

Following Luttinger || Eq. (p]) can be rewritten as 
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where ( = uj + i5/h and 
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For simplicity, as what follows, we denote (T^iC) by °"^( L<J ) EH the calculation of which we 
will concentrate on keeping in mind that we have a finite imaginary part of the denominator 
in Eq. ([|). Consider a pair of eigenvalues e n and e m . For a suitable contour C in the complex 
energy plane (see Fig. 1) the residue theorem implies 

. d M — = _ 2 „ IM — + (4) 
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where the Zk = €f + i{2k — 1)St (cf is the Fermi energy, ks the Boltzmann constant, T the 
temperature and St = T^ksT) are the so-called Matsubara-poles. In Eq. (|) it was supposed 



that Ni and N 2 Matsubara-poles in the upper and lower semi-plane lie within the contour C, 
respectively, i.e., 

[2Nt - l)*r < 5i < (2iVi + 1)<5 T , (5) 
{2N 2 - 1)5 T <5 2 < {2N 2 + 1)5 T , (6) 

Eq. (H) can be rearranged as 

- /M -U ^. .w /W .. + (7) 
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Similarly, by choosing a contour C as shown in Fig. 2 the following expression, 

_ . f{e m ) = — l dz ^ + (8) 
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can be derived. Inserting Eqs. (0) and (§) into Eq. (|3|) and by closing the contours at oo and 
-oo, a^ v {u) is given by 
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It is now straightforward to rewrite Eq. (P) in terms of the resolvent IT 
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such that 
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where Tr denotes the trace of an operator. By using the below quantity, originally introduced 
by Butler [|, 

= ~^F Tr [J il G{z l )rG{z 2 )\ , (12) 
for which the following symmetry relations apply, 

<7vn{z2, Z X ) = ^nu{Zx, Z 2 ), ^(^, Z* 2 ) = Uv^Zx, Z 2 )* = O ^ u (z 2 , Zx)* , (13) 

<Jh U (lj) can be written as 
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which because of the reflection symmetry for the contours C and C (see Figs. 1 and 2) and the 
relations in Eq. ( |i~3"D can be transformed to 



j> dzf{z)a llv {z + hu + z<5, z) — (j> dzf(z)a^ u (z — hu + iS, z)\ (15) 
m 

-2iS T {a^ u (z k + hu + i5,z k ) + a^ u (z k - hu + i5,z k )*} . 

k=-N 2 +l 

Eq. (|15D displays the central result of the present work: it compactly expresses the optical 
conductivity tensor in terms of contributions of a contour integral and those due to Matsubara 
poles. The occurring quantities of the type a^zx, z 2 ) are now exactly of the form such that 
substitutional disorder (CPA) for layered systems |13|, Q can be treated. It should be noted 
that in the case of site-diagonal terms, see [TJ]], in principle for the evaluation of the contour 
integrals also the 'irregular' part of the Green's function has to be included. Since the inho- 
mogeneous CPA for semi-infinite systems and the corresponding Kubo-Greenwood approach to 



electric transport is discussed at length in [14|, no further discussion with respect to layered 
systems is needed. 



3 Integration along the real axis: the limit of zero life- 
time broadening 

In this Section we give the relationship of Eq. ([15]) to formulations existing in the literature. 
For this reason the contour C is deformed to the real axis such that the contributions from the 
Matsubara poles vanish. By using the relations in Eq. flilf ) Eq. (p~5|) reduces to 



/oo 
def(e) [cf^(e + hu + i5, e + zO) — a^(e + hu + i8, e—iO)] (16) 
-oo 
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where we introduced up- and down-side limits for the resolvents [12[]. By taking the limit 5^0 
Eq. ([TJ]) becomes equivalent to Eq. (5.15) of Ref. || for q=0, 
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Note that by shifting in the second term of Eq. fll8|) the argument of integration by hut the 
hermitean part of a^ v {uj) can be expressed as 

1 1 f°° 

CiM") = + <r»M*) = -4^7 y_ oo de (/( e ) - /( e + M) x (19) 

Tr {j"{G + (e + ^)-G~(e + twj)} J u {G + (e) -G~{e)}} . 
Since quite clearly <7i iAti ,(0) = 0, from Eq. (Q) we get Tii ;lxv (u) = ax )liV (u)/u as used in Ref. [ p~Qf| . 

In order to obtain the correct zero frequency formula one has to take first the ui — > limit of 
Eq. (jlBf) and then, after inserting into Eq. @, the 5 — > limit has to be performed. Making 
use of the analyticity of the Green's functions in the upper and lower complex semi-planes this 
leads to || 

£,,(0) = -^ /^de/( e )Tr|^^^^[G + ( e )-G-( e )]-^[G + (e)-G-( e )]^^^| , 

(20) 

which can be integrated by parts. The corresponding expression at T = for the diagonal 
elements yields Butler's original formula 0, 

S w (0) (r=0) = -J-Tr{j»[G+(e F )-G-(e F )}J»[G+(e F )-G-(e F )}} . (21) 



4 Practical evaluation of the contour integrals 

Returning to Eq. (|T5|) we are left with the task to evaluate the occurring contour integrals in 
a manner suitable for computational purposes. This is in particular demanding because, in 
principle, integrations from —00 to 00 are involved. Recalling Fig. 1, it is essential to note that 
the integration pathes can be chosen arbitrarily with the only constraint that the one in the 
lower complex semi-plane should lie in the range > Im(z) > —5. It is trivial to show that, 
if the contributions from the core-states are not considered, the contour C can be closed at 
(bottom of the valence band) which lies energetically just below the regime of valence states. 
Consequently, the integrals in Eq. (|15|) can be split up as follows 

jf dzf(z)a^ v (z ±hu + i5, z) = (22) 

roo+iSi poo — 1S2 

/ dzf(z)a /J/ll (z ± Ttu + iS, z) — \ dzf(z)a lxl/ (z±fnu + i5,z). 

Je b +i0 Je b -i0 

Because of the fast decay of the Fermi function f(z) for Re(z) > e F the integrals on the rhs of 
Eq. (p2]) can be terminated at 

£ u = £F + n kBT (n ~ 5 — 10). 



Finally, after manipulating the contributions from the Matsubara poles, see Eq. (0), the optical 
conductivity (or rather the zero wave-number current-current correlation function |IT||) can be 
written as 



<vM = °IM + °IM + , (23) 

where 
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and 
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+Z»v(eF -hu + i[S + (2k - l)S T ],e F + i(2k - 1)S T )* 
+°nv( e F + + i[5 - (2k - l)6 T ],e F - i(2k - l)8 T ) 
+a^(e F -hu + i[5 - (2k - 1)5t], e F - i(2k - l)5 T )*} . 

Preliminary calculations show that by choosing 5 ~ 0.05 Ryd, 5\ ~ 0.1 — 0.15 Ryd and by 
using Gaussian-quadrature a total of 30 - 100 energy points (depending on Uu) is sufficient to 
sample the above energy integrals. Furthermore, since most of these energy points lie either 
below the valence band or far enough away from the real axis (the energy point closest to the 
real axis is at z\ = e F + i5 F , where for T = 300 K 5 F — Q mRyd), the occurring Brillouin-zone 



integrals (see Ref. (Tj]) can be evaluated by using a sufficiently small set of k-points in the 



surface Brillouin zone, thus ensuring excellent numerical stability. 



5 Summary 

In the present paper, by using a finite life-time broadening (or a finite adiabatic switch-on 
parameter) we reformulated the expression for the (linear response) optical conductivity tensor 
in terms of contour integrations such that Green's functions (resolvents) naturally appear and 
finite temperatures are represented by the Fermi-Dirac distribution function. For zero life-time 
broadening the traditional expressions of the conductivity tensor in terms of the Green's func- 
tion are recovered. Finally, a scheme for calculating the contour integrals that occur, in a nu- 
merically efficient way, is suggested. Although primarily applications to ab-initio calculations of 
magneto-optical properties of ordered and disordered thin films in terms of the fully-relativistic 



spin-polarized Screened Korringa-Kohn-Rostoker method |L4], [15[] are currently under way, it 
should be noted that the present method can also be used to calculate other linear response 
functions, just as readily as it can easily be extended to the evaluation of non-linear response 
functions. 
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Figure 1: Contour in the complex plane corresponding to the integration in Eq. (H). 
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Figure 2: Contour in the complex plane corresponding to the integration in Eq. 



